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ABSTRACT 



We present XMM-Newton EPIC observations of the two nearby starburst merger 
galaxies NGC 3256 & NGC 3310. The broad-band (0.3-10 keV) integrated X-ray emis- 
sion from both galaxies shows evidence of multi-phase thermal plasmas plus an un- 
derlying hard non-thermal power-law continuum. NGC 3256 is well-fit with a model 
comprising two MEKAL components (kT— 0.6/ '0.9 keV) plus a hard power-law (L=2), 
while NGC 3310 has cooler MEKAL components (/cT=0.3/0.6keV) and a harder 
power-law tail (T=1.8). Chandra observations of these galaxies both reveal the pres- 
ence of numerous discrete sources embedded in the diffuse emission, which dominate 
the emission above ^2 keV and are likely to be the source of the power-law emission. 
The thermal components show a trend of increasing absorption with higher tempera- 
ture, suggesting that the hottest plasmas arise from supernova-heated gas within the 
disks of the galaxies, while the cooler components arise from outflowing galactic winds 
interacting with the ambient interstellar medium (ISM). We find no strong evidence 
for an active galactic nucleus (AGN) in either galaxy. 
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1 INTRODUCTION 

By virtue of their greatly enhanced star formation rates, 
starburst galaxies (SBGs) are very powerful sources of X- 
ray emission, and consequently exhibit complex spectral 
characteristics reflecting the variety of energetic phenom- 
ena related to the end points of stellar evolution. Starburst 
episodes are transitory events in the lifetimes of galaxies, 
with enhanced SFRs that are unsustainabl e over periods 
greater than ~ 10 s yr (e.g. lHeckmanlll99^ . Galaxies go- 
ing through starburst phases appear to be a common phe- 
nome na in the early Universe (e.g. ISmail. Ivison. fc Blainl 
Il997f) . and local SBGs show characteristics similar to high- 
redshift examples. Therefore, our understanding of the X- 
ray emission processes in local systems is of fundamental 
importance. For example, SBGs are believed to be the prime 
source of "stellar feedback" , whereby supernova- driven out- 
flows (winds) enrich both the galaxy's int erstellar medium 
(ISM ) and the intergalactic medium (IGMI llStrickland et alJ 
2000). The similarity of local and high-redshift examples 
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is further highlighted by the detection of large-scale galac- 
tic winds in UV o bservations of star - forming Lyman bre ak 
galaxies at z ~ 3 dPettini et alJl2000t IPettini et aljl200ll) . 

Mergers of galaxies appear to trigger the most ex- 
treme starburst episodes. Merger systems tend to have ex- 
tremely high infrared (IR) luminosities (luminous IR galax- 
ies (LIRG): Lir > 10 11 Lq and ultraluminous IR gal axies 
(ULIRG): Lir > 10 12 L Q , e.g. lSanders fc Mirabe]|ll996lh the 
brightest of which can even overlap th e bolometric lumi- 
nosities of quas ars and Seyfert galaxies llLipari et alJl20o"ol : 
iLira et al]l2002l and references therein). Many luminous IR 
galaxies are detected in deep-IR surveys, showing that they 
are an import ant constituent of the h igh-redshift (z > 1) 
Universe fe.g. iFranceschini et a l.||200lf) . In the hierarchical 
merger hypothesis fe.g. l r Romrelfi977lh collisions and merg- 
ers of galaxies are thought to be one of the most domi- 
nant evolutionary mechanisms for galaxies, whereby isolated 
disk galaxies merge through close tida l encounters to even- 
tually form elliptical galaxies |K ormend v fc Sanders! Il992l : 
iRead fc Ponmanl I199& iGenzel et al.ll200ll) . X-rav observa- 
tions of nearby systems with obvious signatures of merg- 
ing at other wavelengths provide a window through which 
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to search for the high-energy signatures of their primary 
sources of power i.e. starburst activity or the presence of an 
active galactic nucleus (AGN). 

Prior to the XMM-Newton and Chandra era, X-ray mis- 
sion such as Einstein, ROSAT and ASCA revealed the sig- 
natures of the various processes present in the closest lo- 
cal SBGs. For example, ROSAT X-ray studies of the soft 
(0.1-2.4 keV) emission in the two archetypal SBGs M82 and 
NGC 253 showed that both had substantial extended dif- 
fuse emission from the disks and halo regions that were 
well modelled with thermal plasmas with kT < 1 keV 
JStrickland. Ponman. fc Stevenjri997t iPietsch et a~l]E)0cl) . 
whereas broader band (0.5-10 keV) ASCA studies of these 
and other SBGs showed that their spectra could be de- 
scribed by models comprising one- or two-temperature ther- 
mal plasma components (kT < 1 keV) plus a harder compo- 
nent fitted with either high-temperature plasma or a non- 
thermal power-law (e.g., M82 fc NGC 253.lPtak et alJll997t 
NGC 3310 & N GC 3690. IZezas. Georgantopoulos fc Wardl 
119981: NGC 3256. lMoran. Lehnert fc Helfandlll999l) . 

However, subsequent to the launch of the XMM- 
Newton and Chandra X-ray observatories, we are now able 
to utilize their excellent complimentary spectral and imag- 
ing capabilites to study the detailed spatial and spectral 
characteristics of these X-ray emission processes. The sub- 
arcsecond spatial resolution of Chandra can resolve indi- 
vidual compact X-ray sources and accurately determine 
their positions, whereas the high throughput of XMM- 
Newton can yield detailed spectral information of bright 
point sources and low-surface-brightness extended emis- 
sion. For example, XMM-Newton observations have given 
an unprecedented wide-field view of the spectral proper- 
ties of the di ffuse emission in th e nucleus, disk and halo 
in NGC 253 iPietsch et alJl2001T) . while Chandra observa- 
tions gave detailed images of point sources and diffuse emis- 
sion in_tlM_centrjl_^egk2nsof the galaxy dStrickland et all 
l200d: IStrickland et all 120021). Similar work has also been 
done with M82 (Stevens. Read, fc Bravo-Guerrerol 120031: 
iMatsumoto et alJl200lF ! 

The main X-ray spect ral components of S BGs have re- 
cently been quantified bv IPersic fc Rephaelil <l2002l) . based 
on the observed spectral properties of M82 & NGC 253 
as derived from BeppoSAX observations and a stellar- 
population evolutionary model. They assessed contributions 
from the gaseous and stellar X-ray processes expected in 
SBGs, i.e. from X-ray binaries (XRBs), supernova remants 
(SNRs), Compton scattering of ambient far-infrared (FIR) 
photons off supernova-accelerated relativistic electrons, dif- 
fuse thermal plasma and a compact nucleus in the form of 
a starburst or AGN. Their results show that at low ener- 
gies (< 2keV), the dominant component is expected to be 
a low temperature (kT ^ 1 keV) diffuse plasma resulting in 
part from the galactic wind itself, but mainly from shock 
heating via the interaction of t he hot low-density wind wit h 
the ambient high-density ISM JStrickland fc Stevensll200Cl) . 
At higher energies (2-10 keV), they predict that the spec- 
trum will be dominated by emission from bright neutron star 
XRBs (high- and low-mass secondaries) in the form of a cut- 
off power-law with a possible contributio n from non-thermal 
Comp ton emission or an AGN if present iPersic fc Rephaelil 

However, recent studies of SBGs with XMM- 



Newton and Chandra have shown that, in the absence of 
a bright AGN, the X-ray output of starforming galax- 
ies are dominated (particularly in the hard 2-10 keV 
band) by s mall numbers of bright extra-nuclea r point 
sources (e.g. lKilgard et alJl2002llColbert et alJl2004) . which 
show increasing evidence of being XRB systems, pos- 
sessi ng black hole rather than neutron star prim aries 
(e.s. iFabbiano fc Whitdliool : iMiller fc Colberdl2004l) . The 
brightest of these discrete X-ray sources are the ultralu- 
minous X-ray sources (ULXs), defined as those possessing 
X-ray luminosities > 10 39 erg s _1 , which exceeds the Ed- 
dington limit for accretion onto a IAMq neutron star (see 
section IB~2*1 for further discussion). 

In this paper, we present XMM-Newton X-ray observa- 
tions of the two nearby starburst merger systems NGC 3256 
& NGC 3310, together with complementary published Chan- 
dra data for NGC 3256 and new Chandra results for 
NGC 3310 (Zezas et al. in preparation) to aid our inter- 
pretation of the spectral information from XMM-Newton. 
This paper is set out as follows. In sections 12.11 and 12.21 
we give a summary of the known properties of each galaxy 
based on multiwavelength studies. In section [3] we outline 
the XMM-Newton data reduction methods used, then in sec- 
tion [I] we present the detailed results of the spectral fitting 
procedures and compare these with previous results from 
ROSAT, ASCA and Chandra observations. We discuss our 
results in section in the context of other recent results 
from studies with XMM-Newton and Chandra, and discuss 
how th ese compare as a whole to the starburst spectral tem- 
plate of lPersic fc Rephaelil i2002fl . Our conclusions are sum- 
marised in section HJ Th e distances we adopt in this paper 
are taken from lSanders et al] i2003f) . which assume a value 
of #o=75 kms-^pc" 1 . 



2 THE GALAXIES 
2.1 NGC 3256 

NGC 3256 is a starburst merger system located at a dis- 
tance o f D=35.4Mpc. It is IR-bright (L FIR = 2.7 x 
10 11 L^. |Sanders et al]l2003l) . although it does not reach the 
ULIRG category. Its starb urst natu re has been confir med 
at UV llKinnev et allll995ft. optical jLfoari et al]l2000T). IR 
jGraham et alJll984UDovon et al.lll994l:lMoorwood fc Olival 
Il994h and radio wavelengths jNorris fc Forbes! Il995t) . It 
shows a highly disturbed structure with prominent extended 
tidal tails spanning ~80 kpc, which are believed to be the sig- 
nature of a merger between two gas-r i ch galaxies of roughl y 
equal size jToomre fc Toomr 4 Il97l lEnglish et all |2003» . 
This merger is th ought to be relatively young (5 x 10 s yr, 
iLi'pari et alll2000ft . because within the highly disturbed cen- 
tral region two distinct nuclei have been resolved at ra- 



dio dNorris fc Forbes) 1995). IR jKotilainen et alJll996l) and 
optiral llLrpari^t^lT 200C ) wavelengths. The northern nu- 



cleus has been identified as a pu 
erful ouflowing galactic wind tLfpari et al.ll200oT) . Despite 
previous speculation based on near-IR and radio observa- 
tions that the southern nucleus may host an AGN (e.g. 
iKotilainen et alJl996tlNorris fc Forbesll995h . this seems in- 
creasingly unlikely as no unambiguous evidence for AGN 
activity was found in the rece nt Chandra X-ra y observa- 
tion which resolves both nuclei jLira et alJ |2002fl. although 
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Table 1. Details of the XMM-Newton observations. 



Target 


Observation ID 


Date 


Useful exposure (s) 


Filter 


Source 


count rate (ct s 1 , 0.3-10 keV) 






(yyyy-mm-dd) 


PN MOS 




PN 


MOS 


NGC 3256 


0112810201 


2001-12-15 


6862 11370 


Thin 


0.55 


0.16 


NGC 3310 


0112810301 


2001-05-11 


6728 8764 


Thin 


0.67 


0.19 



a recent radio study bv lNeff. Ulvestad. fc CampiorJ (2003) 
detected strong radio emission from both nuclear sources 
as well an ULX and show that both nuclei have proper- 
ties consistent with radio-loud LLAGN. The age of the star- 
burst component in this galaxy is ~ 5 — 25 x 10 6 yr, which 
mea ns that it must ha ve started long after the merger be- 
gan iLfpari et alJl20o'oT) . This galaxy was also observed with 
earlier X-ray missions, i.e. ROSAT feoller et al.lll992l) and 
ASCA llMoran et al. l ll999l) . and has been found to have a 
high X-ray luminosit y (observed Lx ~ 2 x 10 41 erg s _1 in 
the 0.5-10 keV band. lMoran et al.lll999l) . 

2.2 NGC 3310 

NGC 3310 is a closer but smaller system than NGC 3256, 
located at a distan ce of ~ 19.8 Mpc and cl assified as 
type SAB(r)bc pec ide Vaucouleurs et alj ll99"lT) . It has a 
distu rbed structure, and the starburst (aged ~ 10 7 — 10 s 
yrs, lElmeereen et all 120021 and references therein) is 
thought to have been triggered by a collision between 
NGC 3310 and a dwarf c ompanion during the last ~ 10 s yr 
jBalick fc Heckmanlll98lh . Although no double nucleus has 
been resolved in NGC 3310, further evidence of a merger 
comes from kinematic peculiarities in that the centre of ro- 
tation of the emission line gas does not c oincide with its nu- 
cleus (Ivan der Kruit fc de Bruvnl Il976l : iBalick fc Heckmanl 
ll98lT) . One of the most striking features of the optical 
emission is a prominent "bow & arrow" structure extending 
100 arcseconds northwest of the galax y's centre, which is 
believ ed to be a remnant of the merger iBalick fc Heckmanl 
I l98ll). Its FIR lu minosity (L FIR = 2.8 x 10 10 L Q , 
San ders et alJ [2003) is comparable to that of M82, al- 
though a factor of ten less luminous than NGC 3256. 
NGC 3310 has been studied in nu merous multiwave- 
lengt h o bservations, inclu di ng opt ical JBahcl^^^Ieckmai] 
19811: IPastoriza et alJ Il993l iMulder fc van Drie 
1996), ultraviolet (UV) J van der Kruit fc de Bruvnl 
197tj: iMeurer et all 11994 ISmith et alJ Il99fj) . IR 
|Telesco fc Gatlevl Il984t IPastoriza et all Il993t> . Ha 
iMulder fc van Drieli Il99& IConselice et all l200dl and 
radio jBalick fc Heckmanl ll98lF1 vmi^erKniit^^^^^ru^J 
| 1976[: IMulder. van Driel. fc Braind Il995t iKreeel fc Sancisi 
2001). Early optical and Ha images showed that the bright 
inne r regions are dominated b y an open spiral arm pattern 
(e.g. IBalick fc Heckmanl Il98lh . and the inner part of this 
region connects to a ~ 30 arcsecond (~ 2.9 kpc) diameter 
starburs t ring, which in turn surrounds a blue compact 
nucleus llKregel fc Sancisll200lf) . Optical and near-IR spec- 
trophotometry observations have shown that the nucleus 
has solar abundances, whereas the c ircumnuclear and dis k 
HII regions show a lower metallicity iPastoriza et al.l|l993l) . 
There is a "jumbo" HII region ~ 12 arcseconds southwest of 
the nucleus, whose size and Ha luminosity is comparable to 



the larges t extragalactic HII regions known (e.g. NGC 5471 
in M101. IBalick fc Heckmanl Il98ll) . Most recently, high- 
resolution Hubble Space Telescope (HST) optical (WFPC2) 
and IR (NICMOS) observations of NGC 3310 show the 
presence of numerous super star clusters in the innermost 
southern spiral arm and circumnuclear ring whose positions 
correlate with radio and Ha peaks, with ages consistent 
with their formation resulting from the ca nnibalism of a 
dwarf galaxy in the last 10 Myr dElmegreen et all2002T) . It is 
important to note that there is no evidence of AGN activity 
in this system at any wavelength . For example, n o broa d 
Ha line is detected by iHo. Filippenko fc Sargernl il997f) . 
which would be expected from an AGN with the hard 
(2-10 keV) X-ray luminosity measured from ROSAT and 
ASCA observations assum ing the Lx I Lhc re lation of 
lElvis. Soltan fc Keel il984h (c.f. lZezas etaHll998ft . 



3 OBSERVATIONS, DATA REDUCTION & 
ANALYSIS TECHNIQUES 

3.1 XMM-Newton Observations & Reduction 

The details of the XMM-Newton observations are shown in 
Table[U NGC 3256 was observed for 17 ks on the 15th Dec, 
2001; NGC 3310 for 19 ks on the 11th May, 2001. During 
both observations the EPIC MOS-1, MOS-2 & PN cam- 
eras were operated with thin filters in "Prime Full Win- 
dow" mode. The event lists were pipeline-processed using 
the SAS (Science Analysis Software) v5.2.1 (NGC 3256) and 
v5.1.0 (NGC 3310) and all data products (spectra, images 
& lightcurves) were created using SAS v5.4.1 and the EPIC 
calibration as of March 2003. Full-field light curves were ac- 
cumulated for all exposures in the two observations to check 
for high background intervals of soft proton flares. In the 
case of NGC 3256, there were numerous small flares through- 
out the observation and a large section of high level flaring 
for the last ~3ks. Time intervals with count rates > 50 ct 
s" 1 (PN) and 15 ct s _1 (MOS) were cut from subsequent 
data analysis, leaving a net good time for each camera of 
~llks (MOS 1 & MOS 2) and -7ks (PN). For NGC 3310, 
there was a high level of flaring throughout almost the en- 
tire observation. In order to extract any useful information 
from the data, we accepted time intervals with relatively 
high background count rates (< 140 ct s -1 [PN] & < 30 ct 
s _1 [MOS]), providing us with relatively clean data below 
6keV. This left net good times of ~9ks (MOS 1 & MOS 2) 
and ~7ks (PN). 

3.2 XMM-Newton Data Analysis 

The X-ray spectra of NGC 3256 & NGC 3310 were ex- 
tracted using circular regions of radii 45 and 40 arcseconds 
respectively. The background regions were taken close to 
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Figure 1. [Top] XMM-Newton EPIC soft (0.3-2 keV, loft) & hard (2-10 keV, right) X-ray images of NGC 3256. The images are the sum 
of all 3 detectors, and are convolved with a 1 pixel (4 arcsecond) HWHM 2-D Gaussian mask and displayed with a logarithmic scale. 
The intensity minima/maxima correspond to surface brightnesses of 3 X 10~ 5 /3 X 10 — 3 ct s~ 1 (soft) and 3 X 10~ 5 /8 X 10 — 4 ct s _1 
(hard), and the contours correspond to 2 X 10~ 4 / 1 X 10~ 3 ct s _1 (soft) and 8 X 10~ 5 /5 X 10~ 4 ct s" 1 (hard). [Lower left] Chandra (0.3- 
8 keV) linearly scaled image of the (dot-dashed) boxed region shown in the XMM-Newton image, adaptively smoothed with a 3<r local 
significance threshold. [Lower right] Chandra X-ray contours corresponding to surface brightnesses of [0.4,0.7,1.1,1.8,2.1,3.2,4.3] xl0 — 4 
ct s" 1 overlaid on an XMM-Newton OM UVW2 (A centr(l ;=2120A) image of NGC 3256. 



each source, using as large an area as possible with approx- 
imately the same DET-Y distance as the source region (to 
ensure similar low-energy noise subtraction). We used the 
SAS task ESPECGET with standard filtering to simultan- 
iously extract source and background spectra as well as cre- 
ate response matrices (RMFs) and ancillary response files 
(ARFs) for each source. The resulting spectra were binned 
to a minimum of 20 counts per bin in order to optimise the 
data for x 2 statistics. 

The spectral analysis of the XMM-Newton data in this 
study has been performed using XSPEC 11.1/11.3. All errors 
are given at the 90 percent confidence level unless stated 
otherwise. For NGC 3256, the MOS and PN spectra have 
been fitted simultaneously in the 0.3-10 keV band, whereas 
we have only fitted the NGC 3310 data in the 0.3-6 keV 



range as the high-energy data were heavily contaminated 
by soft proton flaring. We have also included a free normal- 
isation constant in the models to account for differences in 
the flux calibration of the three EPIC cameras, which differ 
by < 15 percent in practise. 

Following previous work in this field, the two spectral 
models we use to fit the data are a power-law continuum 
to represent the combined emission from compact sources, 
and an optically-thin thermal plasma (MEKAL) to model 
the thermal components expected in SBGs such as galactic 
winds and SNRs. For this study we have chosen to fix the 
abundance parameter in the MEKAL models to solar val- 
ues, as previous results with ROSAT and ASCA data have 
shown that fitting complicated multi-temperature thermal 
plasmas with this type of simple spectral model tends to re- 



© 2003 RAS, MNRAS 000. UTO 



XMM-Newton observations of NGC 3256 & NGC 3310 5 



Table 2. Spectral fitting results for NGC 3256. 



PL 



N H ° 



Mi 

N H a kT(keV) 



M 2 x 2 /d°f 
Ntf a kT(keV) 



Fx b 



Lx c 

Obs Unabs 



Model 1: M+PL (wabs*po+wabs*mekal) 



q+0.25 



0.61 



+0.02 



Model 2: M+M+PL* (wabs*po+wabs*mekal+wabs*mekal) 



1+0.27 



-+0.04 



Model 3: M+M+PL* (wabs*po+wabs*mekal+wabs*mekal) 



.+1.05 



Gal 0.57 4 



9.47 



+ 1.53 



0.73 1 



0.85 1 



387.1/300 9.71 



+0.45 



355.7/298 9.91 



+0.42 



1.46^ 



1.49" 1 



343.3/297 10.00+°^ 1.50j;o24 



3.78 



+0.10 



5.74 



4.53 1 



+0.06 



Notes: Spectral models: PL=power-law continuum model and M=MEKAL thermal plasma (solar abundances). "Absorption 
column in units of 10 21 cm -2 . b Observed fluxes in the 0.3— lOkeV band, in units of 10 — 13 erg s — 1 cm -2 . c Observed and 
unabsorbed luminosities in the 0.3— lOkeV band, in units of 10 41 erg s — 1 (assuming a distance of 35.4 Mpc). fSame hydrogen column 
as applied to the power-law spectral component. *Models with the cool MEKAL Njj component fixed at the Galactic value. Model 
parameter errors correspond to 90 percent confidence limits for 1 parameter of interest. The best-fit model is shown in bold. 



suit in unre alistically l ow me tal abundances of ~0.05-0.3Z© 
iStrickland fc Steve ns 2000 and references therein) . 



4 RESULTS 



4.1 NGC 3256 



In Figure Q we show XMM-Newton EPIC images of 
NGC 3256 in soft (0.3-2 keV, top left) and hard (2-10 keV, 
top right) bands, plus an adaptively smoothed archival 
Chandra image of the central region [lower left]. No X-ray 
structure is resolved in the main body of the galaxy in the 
XMM-Newton data, but one discrete source is detected to 
the southeast (labelled XMM-1, see section 14.1.21 . How- 
ever, the Chandra images do resolve the bulk of the emis- 
sion into several discrete sources embedded within diffuse 
e mission. The resul ts of the Chandra analysis are presented 
in iLira et alj ll2002l) . where it is shown that the integrated 
emission from the galaxy is best fit with a model comprising 
two MEKAL components plus a power-law continuum. This 
study shows that the hard power-law component with an 
index of F ~ 2 is mainly (75-80 percent) due to the four- 
teen compact discrete sources (including the two galactic 
nuclei), which contribute > 20 percent of the total emis- 
sion of the galaxy in the 0.5-10 keV range. [Lira et all (l2002f) 
speculate that the remainder of the hard emission is likely 
to come from unresolved XRBs and supernovae distributed 
throughout the galaxy, and that one of the thermal compo- 
nents (0.6 keV) is the signature of the hot superwind, while 
the other (0.9 keV) arises from SNR heated gas within a 
much smaller volume near the starburst nucleus. Figure 
[lower right] shows the Chandra X-ray contours overlaid 
on an XMM-Newton Optical Monitor (OM) UV image of 
NGC 3256 of the same scale. This illustrates that the UV 
emission, which typically arises in starforming regions, cor- 
relates well with the Chandra X-ray data, although a sec- 
tion of the UV emission is obscured by the presence a large 



dust lane in the southwest of the galaxy (see Figure 10 in 
ILira et alJliool) . 

Since we do not spatially resolve any X-ray structure 
in the centre of NGC 3256 with XMM, we can only fit the 
spectrum of the integrated emission. We began by fitting 
the data with single- and two-component spectral models 
with combinations of power-law and MEKAL components, 
but we rejected these due to unacceptable values of chi- 
squared (x„ > 2), although a M+PL model with equal 
absorbing columns did yield a reasonable fit as shown in 
Tabled In order to improve this, we went on to fit the data 
with three-component models (M+M+PL) to better repre- 
sent the multi-temperature gas expected in this system. We 
also included separate absorbing Nh columns for each of the 
three components, but fixed the absorption in the coolest 
MEKAL component to t he Galactic value (9.5 xl0 20 cm -2 , 
iDickev fc LockmanllT99ol) as fitting this parameter resulted 
in Nh values below the Galactic value. Following the work 
of lLira et alJ i2002T) . we tried a model where both absorbing 
Nh components were fitted independently, and one where 
they were tied together. The results are shown in Table [5] 
the best-fit model (shown in bold) is the one with inde- 
pendent absorbing columns with fcT=0. 57/0.85 keV, r=1.98 
and xt = 1-16, which is a > 99.9 percent improvement over 
the M+PL model according to the F-test statistic. The PN 
and MOS spectra are shown with the best-fit model in Fig- 
ure |21 and Tabled lists the observed fluxes, unabsorbed lu- 
minosities and overall percentage contributions to the total 
of each from the three spectral components in the XMM- 
Newton fit. 

The errors on the fit parameters in Table [5] corre- 
spond to 90% confidence limits for 1 parameter of inter- 
est (Ax 2 =2.71). However, the temperatures and absorption 
values for multiple thermal components in models such as 
these are strongly correlated. The absorption column as- 
sociated with the cool thermal component is fixed in the 
best-fit model, but more realistic errors for the warm com- 
ponent i.e. for 2 parameters of interest for a 90% con- 
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Table 3. Fluxes, luminosities and percentage of total emission 
from the three components of the best-fit M+M+PL model of 
NGC 3256. 



Component 


Fx 


% of total 


Lx 


% of total 


Cool MEKAL 


2.50 


25 


0.53 


12 


Warm MEKAL 


3.35 


33 


3.16 


70 


Power-Law 


4.15 


42 


0.84 


18 


Total 


10.00 


100 


4.53 


100 



Notes: Observed fluxes in the 0.3— lOkeV band, in units of 
10 — 13 erg s — 1 cm -2 . Unabsorbed luminosities in the 
0.3-10 keV band, in units of 10 41 erg s — 1 (assuming a 
distance of 35.4 Mpc). 



fidence limit (Ax 2 =4.61) are practically unchanged with 
N H = 8.0- 11.4 x 10 21 cm -2 and kT=0. 74-0.95 keV, demon- 
strating that these parameters are well constrained. 

The PN data also show possible line emission at ~6- 
7keV, which could be a result of either neutral 6.4 keV Fe-K 
emission from an AGN or helium- like Fe 6.7 keV emission 
from a population of type lb /Ha SNRs residing in regions 
of star- formation fe.g. lBehar et al.ll200lll . We have therefore 
fitted a narrow gaussian line to the power-law continuum of 
the best-fit model, which results in an unconstrained line 
energy of ~ 6.5 keV with an equivalent width of 250 eV. 
Although the addition of the gaussian component only im- 
proves the fit statistics at the la level, the 90 percent upper 
limit of the line normalisation allows us to constrain an up- 
per limit of 1045 eV for the equivalent width of the line. We 
will discuss the implications of this possible line emission in 
section 



4- 1.1 Comparison with previous work 

Reassuringly, the XMM-Newton results are remarkably sim- 
ila r to the results of the 30 ks Chandra observation reported 
bv lLira et al.l ^2002), where the best-fit model for the inte- 
grated emission from NGC 3256 is also a M+M+PL model 
with fcT=0.58/0.92keV and F=1.99. In addition, the ab- 
sorbing hydrogen columns are consistent between the two 
datasets within the 90 percent confidence limits. However, 
the observed flux of the Chandra model in the 0.5-10 keV 
range (~ 1.2 x 10 -12 erg cm -2 s" 1 ) is ~20 percent higher 
than the equivalent XMM-Newton model in the same energy 
range (~ 1.0 x 10 -12 erg cm -2 s _1 ). In order to try to un- 
derstand this difference, we have compared the relative con- 
tributions from the three model components in the XMM- 
Newton data (Table[3J with those found in the Chandra data 
based on the model para meters and compo nent normalisa- 
tions quoted in Table 7 of iLira et alJ ^2002). This compar- 
ison shows that, while the summed contributions from the 
two thermal components agree to within ~10 percent, there 
is a much lower contribution from the hard power-law com- 
ponent in the XMM-Newton observation (~65 percent of 
the Chandra flux). This is not an aperture affect, as the 
XMM-Newton and Chandra extraction regions are virtually 
the same size (45 and 40 arcsecond radii respectively). It is 
much more likely that the difference in the power-law flux is 
a real variation in the summed flux of the ULX population. If 
the brightest compact X-ray source (the northern nucleus, 
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Figure 2. XMM-Newton spectra and folded best-fit model 
(M+M+PL) of NGC 3256: PN [top], MOS 1 [middle], MOS 2 
[bottom]. The model components are denoted by dashed (power- 
law), dotted (cool MEKAL) & dot-dashed (warm MEKAL) lines. 



see Table 3 in ILira et alJl2002h were to "switch off", this 
would result in a ~21 percent drop in the power-law flux, 
which would explain the majority of this difference. The al- 
ternative is a situation in which there is a coincidental drop 
from several sources, although this is much less likely. 

The resu lts of the fits t o the A SCA data for NGC 3256 
reported by iMoran et al ] Jl999h showed similar results 
and model parameters, with the best-fit model compris- 
ing two Raymond-Smith (RS) plasma components plus a 
hard power-law (kT=0. 29/0.80 keV, r=1.68). Even though 
the relative contributions of the ASCA model components 
differ from both the Chandra and XMM-Newton results, 
the overall observed flux in the 0.5-10 keV range (~ 1.3 x 
10~ 12 erg cm" 2 s" 1 ) agrees with the observed Chandra flux 
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to within ~5 percent iLira et alJl2002h . though it is ~30 
percent higher than the XMM-Newton flux. 

4-1.2 Discrete Source Properties 

The discrete source in the XMM-Newton data (named 
XMM-1) at 10 h 27 m 55.3 s , -43° 54' 46" is coin cident with 
source 14 in the Chandra data ijLira et alJl2002ft . Since this 
source does not have sufficient counts for spectral fitting, we 
have computed its hardness ratios using the PN counts in 
three differe nt energy bands u sing a similar method to that 
described in lLira et alJ (|S)02) in order to compare its prop- 
erties in both observations. We compute log(Q.S- lkeV)/(l- 
2keV)=-0.26 and log (2-7keV)/(l-2 keV)=-0.32, which, in 
terms of the XMM-Newton instrument response, is equiva- 
lent to a photon index of T ~ 2.5 and an absorbing hydrogen 
column of Nh ~ 5 x 10 21 cm -2 . This is softer than source 
14 i n the Chandra d ata, which shows F ~ 2 (see Figure 
4 m iLira et al.lf2002l) . although the absorbing columns are 
similarly high. 

If we assume a power-law photon index of T=2.5 and 
adopt an absorbing column of Nh ~ 5 x 10 21 cm -2 (as mea- 
sured from the both the XMM-Newton and Chandra hard- 
ness ratios), we get an observed flux in the 0.5-10 keV range 
of Fx ~ 2.5 x 10~ 14 erg s" 1 cm -2 and an unabsorbed lu- 
minosity of Lx ~ 7 x 10 39 erg s — 1 , within ~ 20 percent of 
the Chandra luminosity for the same parameters. This lumi- 
nosity places this source into the ULX regime, although this 
value is highly dependent on the large absorption correction 
esti mated from the ha rdness ratios alone. 

ILira et all i2002Tl speculate that the luminous X-ray 
sources in NGC 3256 are likely to be high-mass XRBs 
(HMXBs). They detect no short-term variability in the 
Chandra data expected from an accreting source, although 
variability on short timescales is not generally detected i n 
Chandra observations of ULXs (c.f. iRoberts et al.l l2004fl . 
The position is coincident with an outer spiral arm of the 
galaxy and offset by ~4 arcseconds from an apparent fore- 
ground star, but given that there is no evidence of either 
significant long- or short-term variability, this source could 
be an obs cured background AG N. Using the log(A r )-log(5') 
relation of lHasinger et alJ feOOll) . we calculate that in an an- 
nular region of ~ 1.2 arcminute radius centred on NGC 3256, 
excluding the 0.45 arcminute radius of the galaxy itself, we 
expect to detect 0.02 background AGN at the 2-10 keV flux 
of this source (Fx ~ 1.6xl0~ 14 ergs" 1 cm" 2 ). This suggests 
that the source is likely to be associated with NGC 3256. 

4.2 NGC 3310 

b igure 01 [top] shows the soft (0.3-2 keV, left) and hard (2- 
6keV, right) XMM-Newton EPIC images of NGC 3310. The 
majority of the X-ray emission comes from the bright cen- 
tral region of the galaxy, corresponding to the nucleus, star- 
burst ring, inner disk, and outer spiral arm. There are also 
two symmetrically-positioned discrete sources of X-ray emis- 
sion to the northwest and southeast of the nucleus (dis- 
cussed below), which we denote as X-2 (south) following 
IRoberts fc Warwick! J2000t) . and X-3 (north) (although X-3 
was not included in that paper due to an error in the origi- 
nal data analysis). X-3 was also detected in the lZezas et alJ 



lll998ft ROSAT study (J103843.2+533107), although X-2 
was not. 

A recent ~30ks Chandra observation of NGC 3310 
(Zezas et al. in preparation) has revealed the presence of 24 
discrete sources with X-ray luminosities ranging from 10 37 - 
10 40 erg s" 1 , most of which rather spectacularly trace the 
central starburst ring (Figure|3] [lower left]). There is also ev- 
idence in the Chandra images of extended diffuse emission. 
Figure [3] [lower right] shows the Chandra X-ray contours 
overlaid on an XMM-Newton OM UV image of NGC 3310 
showing again, as in the case of NGC 3256, that the UV and 
X-ray emission are well correlated. 

We have fitted the MOS and PN spectra for the 
main body of the galaxy (excluding the two discrete 
sources) following the same fitting procedures for the XMM- 
Newton spectra as used for NGC 3256, albeit only in the 
0.3-6 keV range due to the high level of background flaring 
in this observation (see section T3. 21 . The results are shown 
in Table [I] with fluxes and luminosities extrapolated to the 
0.3-10 keV band. For this galaxy we have also fitted the Nh 
column associated with the coolest MEKAL component, as 
leaving this parameter free does not result in unrealistic sub- 
Galactic values as was the case for NGC 3256. 

To begin with, we fitted the data using models where all 
Nh columns were left as free parameters. This method pro- 
duced very good fits to the data with both M+PL (x 2 =0.89) 
and M+M+PL (x 2 =0.87) models, with a statistical im- 
provement at the 98.2 percent level for the addition of 
the second thermal component. However, the fitted power- 
law component had very little intrinsic column (3xl0 20 
cm -2 ) compared to the thermal components (5xl0 21 and 
4xl0 22 cm -2 for the cool and warm components respec- 
tively). Chandra observations of starforming galaxies have 
shown that point sources typically have non-negligible ab- 
sorption e.g. the NGC 3256 discrete source s have columns i n 
the range 10 21 to 10 22 cm" 2 (see Table 3. lLira et ai]l2002T) . 
Additionally, if we assume that the point sources are re- 
lated to the same stellar population as the thermal com- 
ponents, we would not expect to see such large differences 
in their associated column densities. These models also im- 
plied that at energies < 0.5 keV, the emission was dominated 
by the power-law component from point sources, but this is 
inconsistent with soft (0.3-2 keV) and hard (2-8 keV) Chan- 
dra images, which show a large contribution from apparently 
diffuse emission at soft energies. 

In order to achieve a more physically realistic fit, we 
re-fitted the data with models where the power-law ab- 
sorption component was fixed. Firstly, we fixed it to the 
value obtained from the simple power-law plus absorption 
model (A r ^ = 1.19x 10 21 cm" 2 ), and fitted the remainder of 
the spectrum with one and two thermal components (mod- 
els 1 & 2 in TableQJ . Both models gave good fits to the data, 
with a > 99 percent improvement with the addition of the 
second thermal component. The absorption values obtained 
with these fits are much more realistic than our initial fits, 
with Nh values in the range of 0.4-3xl0 21 cm" 2 . 

Secondly we tried a model in which the power-law ab- 
sorption comprised one Nh component fixed at 1 x 10 21 
cm" 2 , and a partial covering absorption component fixed 
at lxlO 22 cm" 2 with a covering fraction of 0.5, i.e. ef- 
fectively some power-law emission experiencing moderate 
columns and some experiencing higher absorption. This re- 
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Figure 3. [Top] XMM-Newton EPIC soft (0.3-2 keV, left) & hard (2-6 keV, right) X-ray images of NGC 3310. The images are the 
sum of all 3 detectors, and are convolved with a 1 pixel (4 arcsecond) HWHM 2-D Gaussian mask and displayed with a logarithmic 
scale. The intensity minima/maxima correspond to surface brightnesses of 5 X 10~ 5 /4 X 10~ 3 ct s — 1 (soft) and 9 X 10 _5 /8 X 10 — 4 
ct s _1 (hard), and the contours correspond to 4 X 10~ 4 / 1 X 10~ 3 ct s _1 (soft) and 3 X 10~ 4 /6 X 10~ 4 ct s _1 (hard). [Lower left] 
Chandra (0.3-8 keV) linearly scaled image of the (dot-dashed) boxed region shown in the XMM-Newton image, adaptively smoothed 
with a 3<r local significance threshold. [Lower right] Chandra X-ray contours corresponding to surface brightnesses of [2,4,6,16] X 10 -5 ct 
s" 1 overlaid on an XMM-Newton OM UVW1 (A centrai =2910A) image of NGC 3310. 



suited in improved fits (models 3 & 4 in Table^J. Again, the 
absorption values measured with these fits are reasonable 
with N H ~ 6 x 10 20 for the cool MEKAL component and 
N H ~ 2 x 10 21 for the warm MEKAL component. We there- 
fore adopt this M+M+PL model (model 4, shown in bold in 
Table|3J as the best-fit model for this data. In this case, if we 
consider more realistic model errors i.e. 90% confidence lim- 
its for 2 parameters of interest (Ax 2 =4.61), the absorption 
and temperature parameters for the thermal components are 
slightly less well constrained than those of NGC 3256. While 
we can only gain an upper limit to the Nh component tied 
to the warm thermal component (Nh < 4 x 10 21 cm -2 ), 
the constraints on the temperature remain unchanged. The 
parameters for the cool component are only slightly less well 



constrained with Nh 
0.30 keV. 



: 0.3-3.3 x 10 21 cm" 2 and fcT=0.14- 



The PN and MOS spectra are plotted with the best-fit 
model in Figure[I] and TableQ>]lists the observed fluxes, un- 
absorbed luminosities and overall percentage contributions 
to the total of each from the three spectral components. 
The power-law continuum is the dominant component over 
the 0.3-10 keV range, contributing ~80 percent of the total 
observed flux and unabsorbed luminosity of the galaxy. We 
note that we were unable to test for the presence of Fe-K 
lines after excluding data above 6 keV. 
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Table 4. Spectral fitting results for NGC 3310. 



PL 

N H a r 



Mi 

N H a kT(keV) 



M 2 x 2 /d°f 
N H a kT(keV) 



Fx' 



Lx c 

Obs Unabs 



Model 1: M+PL (wabs*po+wabs*mekal) 
1.19+ 



1.70 



-0.07 
-0.08 



2.97~* 



0.26 



+0.03 



-0.78 v -^"-0.03 
Model 2: M+M+PL (wabs*po+wabs*mekal+wabs*mekal) 



1.19+ 



1.56 



-o.io 

-0.07 



0.42 



+4.10 



0.24 



+0.04 



-0.39 v -"-0.( 
Model 3: M+PL (wabs*pcfabs*po+wabs*mekal) 



0.86 



+2.00 
0.45 



1/10* 



1.98 



-0.07 
-0.07 



4.17 



+0.74 
0.83 



0.24 



+0.03 

0.02 



0.63 



+0.07 
0.05 



Model 4: M+M+PL (wabs*pcfabs*po+wabs*mekal+wabs*mekal) 



1/10* 1.84t°;Jl 



0.61 



+ 1.84 
-0.13 



0.25 



+0.04 
0.08 



2.19 



-2.08 

-1.06 



0.62 



-0.06 
-0.05 



322.8/348 



309.6/348 



17.16 



+0.67 
1.84 



288.2/345 18.13 



+0.30 
2.05 



16.53 



-0.70 
-2.07 



283.1/345 17.30t?;|| 



0.81 



r0.03 
-0.09 



0.85 



+0.01 
0.10 



0.78 



0.81 



-0.03 
-0.10 



+0.02 

0.09 



1.40 



+0.06 
0.47 



0.99 



+0.01 
0.09 



2.21 



-0.11 
-0.89 



1.26" 1 



Notes: Spectral models: PL=power-law continuum model and M=MEKAL thermal plasma (solar abundances). "Absorption 
column in units of 10 21 cm~ 2 . b Observed fluxes in the 0.3-10 keV band, in units of 10~ 13 erg s _1 cm -2 . c Observed and 
unabsorbed luminosities in the 0.3— lOkeV band, in units of 10 41 erg s" 1 (assuming a distance of 19.8 Mpc). ' Nh component fixed. 
*Nh components fixed at 10 21 cm -2 and 10 22 cm -2 for wabs and pcfabs components respectively. Errors correspond to 90 percent 
confidence limits for 1 parameter of interest. The best-fit model is shown in bold. 



Table 5. Fluxes, luminosities and percentage of total emission 
from the three components of the best-fit M+M+PL model of 
NGC 3310. 



Component 


Fx 


% of total 


Lx 


% of total 


Cool MEKAL 


1.50 


9 


1.09 


9 


Warm MEKAL 


2.00 


11 


1.98 


16 


Power-Law 


13.80 


80 


9.48 


75 


Total 


17.30 


100 


12.55 


100 



Notes: Observed fluxes in the 0.3— lOkeV band, in units of 
10 — 13 erg s — 1 cm -2 . Unabsorbed luminosities in the 
0.3-10 keV band, in units of 10 40 erg s" 1 (assuming a 
distance of 19.8 Mpc). 

4-2.1 Comparison with previous work 

In the RO SAT and A S CA o bservations of NGC 3310 re- 
ported by Ze zas et alJ <|l99jf ). the combined spectrum of 
the integrated emission from NGC 3310 was well-fit with 
either a double RS thermal plasma (kT=0. 80/14.98 keV) 
or one comprised of a power-law plus RS thermal plasma 
(r=1.44, fcT=0.81keV), with both components affected by 
a single absorbing column. These parameters are partly 
consistent with the XMM-Newton results; the power-law 
slope is slightly harder and the thermal plasma component 
similar (to within the 90 percent confidence limits) to the 
warm MEKAL component in our best-fit M+M+PL model. 
However, the main difference is that whereas they fitted 
a hot thermal plasma plus a low absorbing hydrogen col- 
umn (Nh ~ 0.17 x 10 21 cm -2 ), the improved quality of 
the data in the XMM-Newton observation at the soft en- 
ergies means that we are able instead to model the soft 
emission with a cool thermal plasma with a higher absorb- 
ing column. There is also a marked difference in the mea- 



sured fluxes and luminosities in ROSAT /ASCA and XMM- 
Newton datasets. Even if we consider the nearly identical 
two-component PL+M/PL+RS models, there is a ~ 40 per- 
cent reduction in the observed flux in the XMM-Newton data 
compared with the ROSAT /ASCA analysis (~3x 10~ 12 
erg s _1 cm" 2 ). It is likely that this is due to the difference 
in size of the extraction regions used. The ROSAT spectra 
were extracted using a 1.5 arcminute radius region and the 
ASCA regions were larger at 2.7 and 5.5 arcminu tes radius 
for the SIS and CIS respectively iZezas et alll998f) . whereas 
we have used a 40 arcsecond radius circular source region to 
avoid contamination from the two discrete sources. 

In order to make a direct comparison of the 
observed fluxes, we have attempted to replicate the 
ASCA /ROSAT results by extracting XMM-Newton spec- 
tra in an aperture with the same radius as the largest 
ASCA aperture (5.5 arcminutes) together with as large an 
annular background region as possible given the size of the 
MOS and PN chip arrays. However, this turns out to be 
impractical, as the XMM-Newton spectra are dominated by 
the high background flux (especially at energies > 2keV), 
resulting in low-quality background-subtracted data which 
cannot constrain the spectral shape, and hence cannot con- 
strain the flux from a larger region. 



4-2.2 Discrete Source Properties 

The two discrete sources to the south (X-2) & north (X-3) of 
the main body of the galaxy have insufficient counts for de- 
tailed spectral fitting (~130 and ~110 in each MOS camera 
respectively). However, to estimate their spectral shapes, 
we have co-added and fitted the MOS spectra (X-3 is lo- 
cated on a chip gap in the PN) with simple power-law plus 
Galactic absorption models. Their XMM-Newton positions, 
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Table 6. Discrete Sources in NGC 3310. 
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Figure 4. XMM-Newton spectra and folded best-fit model 
(M+M+PL) of NGC 3310: PN [top], MOS 1 [middle], MOS 2 
[bottom]. The model components are denoted by dashed (power- 
law), dotted (cool MEKAL) & dot-dashed (warm MEKAL) lines. 



power-law photon indices, observed fluxes and unabsorbed 
luminosities are shown in Table |S| The X-ray luminosities of 
these sources are high (7 x 10 39 erg s _1 & 5 x 10 39 erg s _1 , 
0.3-10 keV), putting them well into the ULX regime. Al- 
though not shown in the OM image ( Figure [3] [lower right]), 
they are both coincident with regions of UV emission, which 
is indicative of star formation activity though there are no 
obvious UV counterparts in this data. Although this XMM- 
Newton observation is very flare-contaminated, we tested 
the variability of each source in the XMM-Newton data by 
deriving short-term light curves. The data from the MOS 
cameras were co-added to improve the signal-to-noise ratio, 
and the data were tailored so that each bin had at least 20 
counts after background subtraction. We performed \ 2 tests 
to search for large amplitude variability against the hypoth- 



Source 


Position (J2000) 




r 


Fx 


Lx 


X-2 
X-3 


10 , '38 m 50.2 a +53° 29' 
10 ,l 38 m 43.3 a +53° 31' 


25" 
00" 


1 4O+0-96 
1 66+ 88 

l.DD_ 75 


1.47 
1.05 


7.00 
5.01 



Notes: Observed fluxes in the 0.3-10 keV band, in units of 10 



erg s 1 cm 

in units of 10 39 erg s — 1 (assuming a distance of 19.8 Mpc) 



Unabsorbed luminosities in the 0.3— lOkeV band, 
erg s" 



esis of a constant count rate, but neither source showed any 
significant variability (> 3a). In order to search for longer- 
term variability, we have made a direct comparison between 
their observed fluxes in this XMM-Newton observation and 
the 1995 ROSAT HRI 41 ks observation, both in the 0.3- 
2.4 keV band as this is the energy range covered by both 
mission. We derived fluxes for the ROSAT HRI observa- 
tion with webPIMMS usin g the observed count rate from 
iRoberts fc Warwick! i2000T> for X-2, and a new count rate 
estimate for X-3 from the archival HRI data, assuming the 
absorbed power-law continuum slope measured in the XMM- 
Newton observation and normalized to the 0.3-2.4 keV band. 
During this ~6 year period, the flux for X-3 has decreased 
by ~57 percent, while that of X-2 has increased by ~65 
percent. Variability on t his timescale is consi stent with that 
found for other ULX fc.f. lRoberts et al.l2004l) . and supports 
the scenario that these sources are single black-hole XRBs 
rather than concentrated regions of star formation, as multi- 
ple sources such as a cluster of XRBs would not be expected 
to all change their output simultaneously. 



5 DISCUSSION 

Our results are consistent with other local SBGs stud- 
ied with XMM-Newton and Chandra, e.g. NGC 253 



JPietsch et alJl200ll: IStrickland et alJl2000]: Istrickland et al 



I 2002F M82, iStevens et all 120031: iMatsumoto et al 
120011) and the closes t merger system NGC 40 38 /4039 
(The Antennae), iFabbiano. Zezas. fc Murravl l200ll: 
i Fabbiano et al J 120031; IZezas fc Fabbianol 120021 IZezas et^l] 
l2002al: IZezas et al.ll2002bl) . The X-rav emission from these 
systems are well-fitted with thermal components ranging 
between ~0.2-0.9keV plus harder (> 2keV) emission 
dominated by power-law continua from discrete sources. 

Even with the relatively short exposures of these XMM- 
Newton observations, especially in the case of NGC 3310 
where there was substantial background flare contamina- 
tion in the data, we are able to improve the definition of the 
plasma temperatures and power-law slopes compared to ear- 
lier ROSAT and ASCA data for these galaxies. Additionally, 
the information from the spatially-resolved Chandra data 
has allowed us to deconvolve this spectral information and 
properly assess the contributions from point sources and dif- 
fuse components. 

An interesting result of this study are the different tem- 
peratures of the thermal components in NGC 3256 and 
NGC 3310. To investigate the reasons for this, we have 
summarized the known properties of each system in Ta- 
blcQ We have also included details for the Antennae system 
(NGC 4038/39), as this is the best example of a close merger 
system of two equal-mass galaxies. For the purposes of this 
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Table 7. Comparison of the properties of NGC 3256, NGC 3310 and the Antennae (NGC 4038/9). 



Parameters 


NGC 3256 


NGC 3310 


Antennae (NGC 4038/9) 


Refs 


Distance (Mpc, Ho=75 km s _1 Mpc) 


35.4 


19.8 


19.3 


1.2 


Merger Type 


Major merger - two 
equal mass galaxies 


Merger with dwarf 
companion 


Major merger - two 
equal mass galaxies 


3,4,6 


Merger Age 


~ 5 x 10 s yr 


< 10 s yr 


2 - 5 X 10 s yr 


3,4,5,6 


Starburst Age 


5 - 25 X 10 6 yr 


10 7 - 10 s yr 


~ 5 - 100 X 10 6 yr 


3,5,7 


Dynamical Mass (1O 1O M ) 


~ 5 


~ 2.2 


~ 8 


3,8 


Lfir {LqY 


1.8 x 10 11 


2.0 x 10 10 


2.9 x 10 10 


1 


SFR /H (Mq yr- 1 )* 


32.6 


3.6 


5.2 




L x * ( erg s _1 ) 


4.5xl0 41 


1.3xl0 41 


6.7xl0 40 




Mi* (kT/keV) 


0.6 (12%) 


0.3 (9%) 


0.4 (17%) 




M 2 * (kT/keV) 


0.9 (70%) 


0.6 (16%) 


0.7 (14%) 




r* 


2.0 (18%) 


1.8 (75%) 


1.6 (69%) 





References: (1) Sanders et al. (20031 : (21 Zezas et al. 2002a: (3) Li'pari et al. (2000') : (41 Balick fc He ckman (1 981') : 

(51 Whitmore et al. (1999 1 : (61 Barnes (1988 - ! : (71 Elmegreen et al. (2002) : (81 Kregel fc Sancisi (20011 . tpiR luminosities calculated 

using FIR fluxes derived from the 60/x and lOOfi IR flux densities of Sanders et al. (2003) using the relation of 

Hclou. Sojfe^^^^gwar^gbinson ^8^: FIR = 1.26 X 10~ 11 (2.58560^ + Sl00^) erg cm -2 s _1 . tSFRs determined using the relation 
of Kennicutt (1998): SFH=Lfir/2.2 X 1O 43 M0 yr -1 . * Unabsorbed luminosities in the 0.3-10 keV band and spectral parameters from 
the best-fit models. Percentage contributions from the components to the total unabsorbed luminosity are shown in parentheses. 



comparison, we selected data from the Chandra archive (Obs 
ID 3042) to provide us with a deep observation (~ 70 ks) un- 
affected by background flaring. We extracted a cumulative 
spectrum from all X-ray sources in a central 1 arcminute 
diameter circular region using the CIAO tool ACISSPEC 
(which provides RMFs and ARFs weighted over the differing 
calibration across the extraction zone), and used adjacent 
source- free regions to provide a background measurement. 
We fitted the spectrum in the 0.3-10 keV range with the 
same spectral models used for NGC 3256 and NGC 3310, 
and the data are well- fit with a power-law (F=1.6) and 
two thermal plasmas (fcT=0.40/0.74 keV), with no absorp- 
tion required ab ove the Galactic value (3.4 x 10 20 cm -2 , 
IStark et alll992ft . The power-law component can be entirely 
account ed for by the system s' population of compact point 
sources l|Zezas et alJl2002all . and the plasma temperatures 
are consistent with the detailed sp atial and spectra l analy sis 
of the diffuse emission reported in lFabbiano et all ll2003fi . 

On inspection, albeit on the basis of only three sys- 
tems, there does appear to be a correlation between the 
type of merger, the age of the starburst component and 
the remaining properties {Lfir, Lx & plasma tempera- 
ture). NGC 3256 is a major merger with a young star- 
burst, and displays high levels of FIR and X-ray emis- 
sion together with hot plasma temperatures. Conversely, 
NGC 3310 is a minor merger with an older starburst com- 
ponent, has lower FIR and X-ray luminosities and dis- 
plays cooler plasma temperatures. The Antennae system, 
on the other hand, has in termediate aged star- formation 
spread over ~ 5 — 100 Myr iWhitmore et aljll999l) . and dis- 



plays plasma temperatures between those of NGC 3256 and 
NGC 3310. However, although the Antennae is an equal- 
mass merger like NGC 3256, its FIR and X-ray luminosities 
are a factor of ~ 6 lower; given that the Antennae has the 
greater mass of the two, this suggests a decrease in the power 
of the starburst with time that is independent of the total 
mass of the system. There is also a difference in the relative 
contributions from the spectral components to the total in- 
trinsic X-ray luminosities of the systems. While NGC 3310 
and the Antennae are dominated by the power-law compo- 
nent, the intrinsic emission in NGC 3256 is dominated by 
the hot thermal plasma due to its high absorption column. 
The ratio of intrinsic X-ray and FIR luminosities is also dif- 
ferent between the the galaxies; while NGC 3256 and the 
Antennae both have Lx / Lfir ~ 6 x 10~ 4 , NGC 3310 has a 
greater relative X-ray luminosity with Lx / Lfir ~ 2 x 10 -3 . 

These properties pose a number of interesting questions. 
For example, does the apparent correlation between plasma 
temperature and starburst age simply indicate that the hot- 
ter phases of the ISM are escaping the systems at earlier 
times after the onset of the starburst episodes? Or could 
the higher masses and hence deeper gravitional potential 
wells of NGC 3256 and the Antennae mean that the hot- 
ter components are retained in these systems for longer? Is 
the highly absorbed hot thermal plasma in NGC 3256 also 
a consequence of the younger age of the starburst i.e. has 
insufficient time elapsed for the neutral absorbing material 
to be blown away from the dense central regions? Larger 
samples of high-quality X-ray data from merging galaxies 
would be required to investigate these trends properly. 
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5.1 Thermal X-ray emission 

The soft emission in SBGs is thought to arise primarily in 
hot gas from both the interaction of supernova-driven galac- 
tic winds with the ambient ISM and from SNRs themselves. 
Thermal components are al so seen in more di s tant exam- 
ples of merger systems, e.g. iFranceschini et all |2003) have 
recently reported the presence of soft thermal plasma emis- 
sion in all galaxies in a survey of ten ULIRGs. 

The soft emission in both NGC 3256 and NGC 3310 is 
well-modelled by multi-temperature thermal plasma mod- 
els. Both galaxies show a trend of increasing absorbing col- 
umn densities with increasing plasma temperature, simi- 
lar to that found in other SBGs. Although we cannot spa- 
tially resolve the differen t emitt ing regions with the XMM- 
Newton data, iLira et all l)2002f) were able to show that in 
the Chandra data, the hot 0.9 keV component dominates in 
the inner ~2 kpc of NGC 3256, while the cooler 0.6 keV com- 
ponent dominates at radii > 2.5 kpc (adjusted to a distance 
of 35.4 Mpc). This evidence implies that the absorbed hot 
phase of the ISM arises from supernova-heated gas possibly 
within the galactic disk, while the cooler component with 
less absorption may arise in outflowing galactic winds. The 
thermal components in NGC 3310 have characteristics con- 
sistent with this picture, although spatially resolved spec- 
troscopy are required to pin down the different emitting re- 
gions. 

As this study demonstrates, particularly in the case 
of NGC 3310 where we have successfully modelled a two- 
temperature thermal plasma for the first time, the increased 
photon statistics and hence spectral quality we can acquire 
with the current generation of X-ray telescopes is reveal- 
ing increasingly complex temperature profiles of the thermal 
emitting gas. However, we must bear in mind that this type 
of modelling with distinct thermal plasma temperatures is 
likely to be over-simplification of the real state of the hot 
gas in these systems; instead a phase continuum of states 
with te mperatures between T ~ 10 s — 10 8 K is more likely 
to exist Strickland fc StevendfeoOOft . 

5.2 Hard X-ray emission 

The Chandra data for these galaxies has resolved the ma- 
jority of the hard power-law flux in both of these galaxies 
into a population of luminous compact sources embedded 
in the diffuse emission (a full analysis of the Chandra data 
will be presented in Zezas et al. in preparation). This agrees 
with other recent studies of SBGs, where large numbers of 
ULXs have been found to be associated with regions of 
curre nt star formation in e.g. The C artwheel llGao et alJ 
l2003h : The Ant ennae (NGC 4038/3 9; IZezas et al]l2002blk 
NGC 4485/90 iRoberts et al] 120021 : Arp299 (NGC 3690; 
IZezas. Ward, fc MTrravll2003f) . suggesting that they are as- 
sociated with young stellar populations and are therefore 
likely to be HMXBs g oing throug h a phase of thermal- 
timescale mass transfer dKing|l2004l) . Such studies serve to 
emphasize the dominant contribution of a small number 
(~10) of probable black hole XRBs to the flux of the power- 
law continuum component in the most luminous SBGs. 
We note that this disag rees with the suggested model of 
IPersic fc Rephaehl {2002), where the power-law continuum 
should originate in a large number of neutron star bina- 



ries (> 500 to reach ~ 10 41 erg s _1 if they are Eddington- 
limited) . Future composite models for SBGs should therefore 
consider this subtle change in the origin of the power-law 
component. 

Although these few sources dominate the emission at 
high energies, there is also likely to be a contribution from 
hotter gas that we cannot de t ect w ith this data. In the 
case of NGC 3256, ILira et all i2002l) have shown that af- 
ter point source subtraction, the hard component is best 
modelled with a ~4keV thermal plasma, though there is 
likely to be a contribution from unresolved point sources. 
Similary in NGC 253, emission lines from thermal plasma 
component > 5 keV are detected in the XMM-Newton RGS 
(Reflection Grating Spectrometer) and a hot plasma of 
~6ke V is required in th e EPIC spectra of the nuclear re- 
gion iPietsch et alj l2001L 

5.3 Starburst/AGN activity 

We find no strong evidence for the presence of an AGN in 
either NGC 3256 or NGC 3310. Hidden AGN activity can 
be detected at X-ray energies by the presence of a highly ab- 
sorbed (Nh > 10 22 cm -2 ) hard component and/or a strong 
neutral Fe-K line at ^6.4 keV originati ng in the nucleus. 
In the XMM-Newton ULIRG survey of IFranceschini et alJ 
(2003), roughly half of the systems studied showed evidence 
of AGN activity, confirming their composite nature but sug- 
gesting that the starburst component is dominant as the 
thermal signatures were detected in all systems. Evidence 
for AGN activity has also recently been found in nearby 
ULIRGs with the d etection of strong neutral Fe-K emission 
lines e.g. NGC 6 240 feoller et all2003tlKomossa et all2003ft 
and NGC 3690 feallo et al.ll2004l) . 

However, the power-law components in the best-fit 
models of both galaxies do not require such high absorb- 
ing columns (NGC 3256, An ~ 1.5 x 10 21 cm -2 ; NGC 3310, 
Nh ~ 1.2 x 10 21 cm -2 ), and have been shown to origi- 
nate in the discrete source populations. There is marginal 
evidence of Fe-K line emission in NGC 3256, with an uncon- 
strained ~6.5keV line energy and a 90 percent upper limit 
of ~lkeV. This is consistent with the equivalent widths ob- 
served in some Com pton-thick type 2 AGN, which ca n be as 
high as >lkeV (e.g. lMolendi. Bianchi. fc Matdl2003h . How- 
ever, hot thermal starburst components can also give rise 
to h elium-like ^6.7 keV emission from type Ib/IIa SNRs 
(e.g. iBehar et all 1200 llL This type of emission is seen in 
M82 iStevens et alJ I 2003T) and the starburst c ompent of 
NGC 4945 dSchurch. Roberts, fc Warwick! 12002) . However, 
this is only a weak upper limit on the equivalent width of 
a la line detection, and as we are unable to distinguish be- 
tween the two types of line emission in these data, the ques- 
tion of the presence of a hidden AGN remains open. 



6 CONCLUSIONS 

In this paper, we have studied the X-ray spectral proper- 
ties of the two nearby starburst merger galaxies NGC 3256 
and NGC 3310. We have used the broad-band (0.3-10 keV) 
spectra from XMM-Newton together with complementary 
high-spatial-resolution imaging data from Chandra to as- 
sess the various stellar and non-stellar contributions to the 
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X-ray emission. We have shown that the emission < 2 keV is 
dominated by thermal emission believed to arise in the inter- 
action of galactic winds with the ISM and SNRs, which can 
be well-modelled by thermal plasma models with representa- 
tive temperatures of ~0.3-0.9 keV. The emission between 2- 
10 keV is dominated by power-law emission from the bright 
ULX population, with the remainder likely to come from 
fainter unresolved point sources and a ho tter plasma. These 
result s are consistent with the models of IPersic fc Rephaelil 
(2002), albeit with the slight modification that ULXs dom- 
inate the high energy emission, and provide yet further ev- 
idence of the complex nature of SBGs. The more thorough 
understanding of local merger and starbursting systems ob- 
tained in studies such as this will certainly assist in our 
interpretation of the very photon-limited X-ray emission we 
currently detect from high-redshift merger systems. 
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